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hole is computed. 
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1. Introduction 



It has been shown by Bizoh, Chmaj and Schmidt that in odd number of spacetime 
dimensions one can perform a consistent cohomogeneity-two symmetry reduction of the 
vacuum Einstein equations which - in contrast to the spherically symmetric reduction - 
admits time dependent solutions. The key idea was to modify the standard spherically 
symmetric ansatz by replacing the round metric on the 5" with the homogeneously 
squashed metric. By squashing we mean here the construction of a metric for 5"" that 
does not admit full symmetry group SO{n + 1). Instead, the metric is required to 
be invariant under transitive action of a specific subgroup of SO{n + 1). This idea, 
technically called the BCS ansatz, has been used for the stability analysis of the 4+1 
dimensional Taub-NUT soliton, also known as the Kaluza-Klein (KK) monopole [2]. 
The monopole is important in M/String theory, as it is super-symmetric and can be 
used as a base for construction of certain BPS states [3], 0]. The KK monopole was 
found to be stable within the ansatz (for small perturbations). 

One might wonder why the KK monopole is stable. In order to address this issue 
we decided to change dimensionality of the model. This way we could figure out which 
features are dimension dependent. We chose to go from the 4+1 to the 6+1 spacetime 
dimensions. There is an additional advantage of this approach connected with the fact 
that the 6+1 dimensional Taub-NUT soliton (TN) is not super-symmetric. Therefore 
we can also see which aspects of relaxation rely on super-symmetry. We found that 
introduced changes do not destroy stability, as the TN soliton is nonlinearly stable 
against small finite-energy perturbations within the BCS ansatz. Our techniques used 
for stability investigation follow (and extend) the ones used in [2], therefore the reader 
is encouraged to read both papers in conjunction. 

Further research revealed that strongly perturbed TN soliton loses stability and 
collapses, forming a black hole. This also occurred in the lower dimensional model. 
This time, however, the black hole solution is not known in closed form. We examined 
it in detail, proving its existence and computing its asymptotic expansion. This allowed 
us to compute its mass and compare it to the mass of its lower dimensional cousin. 



2. Taub-NUT soliton 



The Taub-NUT soliton in 6+1 dimensions is a static regular nontrivial solution of the 
vacuum Einstein equations. It can be explicitly found in [4]: 



'6+1 - -dt + ^p^p^^i^dp + Pip + 2l)dE, + ^^^2/)2 ^ ' 



where / is a free positive parameter. 



a' ={dp-^sm'^as)\ (2) 
dJ:l = d^' + ^ sin^ e {al + (T^ + cos^ ^ al), 
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the latter being the Fubini- Study metric for CP^ and 

0"! = sin ip sin Odcj) + cos ipdO, 
02 = cos ip sin 6d(j) — sin ipdO, 
a^ = dip + cos Odcp 

are standard Maurer-Cartan forms for SU{2). Parameter ranges are as follows 

< p, < < 27r, 0<e <7r, 
0<ip <47r, < ^ < ^, < p < 27r. 

It is worth noticing that this solution lacks asymptotic flatness in a fashion similar 
to its 4+1 dimensional counterpart. It approaches times a circle of asymptotic radius 
a/2/ near infinity. Near p = 0, on the other hand, it is locally fiat, which can be seen 
by introducing a new radial coordinate s = y/2lp. Then the solution is interpolating 
between two different vacua which justifies the name: soliton. 



3. Preliminaries 

The squashed can be constructed as a coset SU{3)/ SU{2) with SU{3) C S0(6) 
being its symmetry group. The result is widely known and can be easily derived using 
the methods of O [6l [71 [8] . The BCS ansatz has the form 

ds^ = -Ae-^'de + A-^dr^ + r^e-^'^a^ + e^^rfS^), (3) 

where cr^ and dTP^ were defined in ([2]) and A, S, 5 are functions of r and t. The function 
5, which measures squashing resulting from lower symmetry of our description of S^, 
serves as a dynamical degree of freedom and breaks Birkhoff 's theorem for that case. It 
has been chosen in such a way that it does not change spheres' volume. 

Substituting the metric ([3]) into the vacuum Einstein equations we get the following 
system of equations 

A A 1 

A = -— + -(24e-2« - 4e-i2^) - Ar{e^'A-^B^ + AS'^), (4a) 
r 5r 

A = -SrABB', (46) 
6' = -Ar{e^^A-^B^ + B'^), (Ac) 

(e^A- ' - {e-'Ar^'B')' + ^e" V(e-2^ - e"'^'') = 0, (Ad) 

where primes and dots denote derivatives with respect to r and t, respectively. Note 
that the equations are scale invariant. The Taub-NUT solution ([T]) expressed in terms 
of ansatz ([3]) takes the form 

.-lOBo 2l\p + Al) 2{p + Al)p dr _ 

' - ip + 2iy ' (p + 2/)2 %^ ' ' 

ri° = 2/V(p + 2/)^(p + 40. (6) 
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It can be proven that every static regular solution of fl4ati4(ip should exhibit Taylor series 
expansion near r = 

5(r) = 6r2 + 0(/), A(r) = l-86V^ + 0(r^), 5(r) = -46V + 0(r^)(7) 

with a/6 being the scale factor (a rigorous proof can be found in appendix A). It is easy 
to verify that ([5]) possesses the same expansion, with 6 = Thus we identify / with 
the scaling freedom and conclude that up to scaling there exists the unique solution (E]). 



4. Linear stability 

The linear stability analysis proceeds along the same lines as in the 4+1 dimensional 
case. Assuming 

A(r,t) = Ao(r) + Ai(r,t), B{r,t) = B^{r) + B^{r,t), (8) 
5(r,t) = 5o(r) + 5i(r,t), 

and keeping only linear terms in perturbations, we get an equation for a single Fourier 
mode (p(r) is defined by equation ([6])) , Bi{p{r),t) = f (p)e''^*: 

Lv{p) = XMp), (9) 

with 

-2 d 3 , d \ 



where V{p) is a long and complicated expression (given in the appendix B). This 
equation cannot be solved analytically for generic A^. However, in order to demonstrate 
linear stability we only need to know the zero mode. This mode can be deduced from 
the scaling freedom. For 6 near 1 we have 

5o(6r) = 5o(r) + r^Bo{r){b - 1) + 0((6 - l)^). (10) 
dr 

We define Vq = r-^BQ^r). Changing r — )■ p we get 

(p + 50p 
4(p2 + 5p/ + 5/2 

with < p, < /. It can be easily verified that this solves (jO]) for = 0. We see 
that the zero mode is manifestly positive for all p in the domain. This implies by the 
standard Sturm-Liouville theory that the operator L has no negative eigenvalues. Thus, 
the TN soliton is linearly stable within our ansatz. Note, that the zero mode is not 
a genuine eigenfunction because it is not square integrable. Rather, it is a resonance 
sitting at the bottom of the continuous spectrum. 
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5. Warped black holes 

Black hole solutions play the key role in dynamics of a strongly perturbed Taub-NUT 
soliton. As far as we know, in contrast to the 4+1 dimensional case, these solutions are 
not known in closed form in 6+1 dimensions. Thus, from the mathematical point of 
view the very existence of such solutions is not obvious and requires a proof. Following 
Breitenhohner, Forgacs and Maison [9] one can prove the local existence (see Appendix 
A). Global existence proof, as well as analysis of an asymptotic behaviour for large r is 
described below. 

Let us begin by explaining what we mean by a black hole solution in our model. 
If we assume that for a given th > ^ Airn) = and substitute the formal power series 
expansion into the static equations, we obtain: 

A(r)~ -(66-2/^-6-12^) (—-I 
5 \rH 

Theorems presented in the Appendix A guarantee that a local solution with the above 
behaviour for small r — exists and is analytic in P,rH and r — r^. By setting (3 = 
one arrives at the ordinary Schwarzschild black hole with B = 0. Numerical simulations 
show that solution with /3 < plays no role in TN soliton evolution. Below, the global 
solution with local behaviour (fT2|) and f3 > will be called the warped black hole 
(WBH). 



5.1. Global properties 

We now focus on establishing the global existence of WBH. Using logarithmic variable 
r = Inr, static version of the equations f l4aMd[) can be cast into the following form: 

AB" + - [B' (246-2^ - 46-12B) _ g ^^-2B _ g-12B)] ^ (^3^) 

A' = -AA (1 + B'^) + - (246-2^ - 46-125^ ^ (^3^) 
S' = -AB'^ (13c) 

Here prime corresponds to differentiation with respect to r. This is an autonomous 
system, which is a consequence of scaling invariance. Note that the first two equations 
are independent of the third one. 

It can be easily verified that WBH solution obeys the following inequalities in a 
right-open range [th,Ti) where Th = Inrn'- 

< fi' < i, 0<A<1, 0< B. (14) 

Our argument proceeds in two steps. In the first step we show that if the local WBH 
solution preserving f[T^ is given in some right-open interval / = [tq,t^), then the 
solution can be safely extended to / = [To,r+] preserving (fT^ . In the second step we 
prove that the first step is equivalent to global existence of the WBH. 



Nonlinear stability of the Taub-NUT soliton in 6+1 dimensions 



6 



Step I 

Let us define an auxiliary function 

and observe, tliat ^ — e > f{B) for r G / and some finite e > dependent on tlie clioice 
of /. Now, if we take B" < and look at fll3al) we get 

B' > f{B), (16) 

for r & I. Therefore, we can safely assume that B" < for B' > ^ — e. This naturally 
leads to a conclusion that B'{t^) < \. Next, we can put 5' = in (II Sal) to obtain 

B" = - > 0, (17) 

which means that i?'(r+) > 0. We can proceed with similar considerations for the 
function A. Setting A = 1 in fll3&P we get 

A = -AB'^ - i (20 + 4e-i2i? _ 246-'^) < 0, (18) 

which means that A{t^) < 1. Finally, setting A = in (11361) we get 

A' = -(2Ae-^^ - 4e~i2B) > g, (19) 
5 

which means that A{r^) > 0. One should not be concerned with setting A = in 
— 4A(1 + B'"^) because B' is finite in I. Finally, from 5' > in / we get that i? > in 
/. Note, that B{t+) < oo because B'{t^) < oo. 



Step II 

Let us assume that a local WBH solution satisfying (fT4|) in some [th, t+) cannot be ex- 
tended to infinity preserving fll4p . Then, there exists > r+ where (I14p is not satisfied. 
By our construction, the solution must be finite and preserve f fT^ for / = [r+,rs). We 
can now use the results of Step I to get a contradiction - the solution is regular and 
preserves ( !T^ in r^. 



We shall now prove that 5 ~ r for large r. To see this, let us first note that 
B{t) — )■ oo as r — )■ oo. Indeed, we already know from f lT^ that B is monotonically 
increasing. Furthermore, if it was to reach a finite limit as r grows, B' would have to 
approach and f{B{r)) would attain a finite positive limit Yrnir-^oo f {.B{t)) = ^ > 0. 
However, if we reverse a sign in ( IT6|) we immediately arrive at the contradiction because 
whenever B' drops below B" becomes positive repelling B' from 0. 

Next, let us observe that 

A/(5(r)) = B'nB) > 0. (20) 

In the neighbourhood of r^^ = Inr^/ we have B' = f{B) which follows from (fT2|) . As we 
depart towards greater r, B' must drop below f{B) because from (ITBi) we know that 
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initially B" = < ^/(5(r)). It becomes evident, by the same argument, that once 
B' drops below / it can never cross it again. Therefore B' is monotonically increasing 
(from reversed ( |T6l) ) and bounded from above (because / is bounded from above and 
B' < f). We get lim^_i.oo B'{t) = ^ > 0, which implies 

i?(r)~Cr + S (21) 

for large r. S is an integration constant that, as we shall see, corresponds to the scale 
invariance. 



5.2. Exact asymptotic behaviour 

We shall now focus on establishing the large r series expansion for WBH solutions. 
First, notice that substituting (12T|) into fll3al) and neglecting e~^^^ yields the consistency 
condition 

(4C - l)e-2^ = (22) 

which fixes C at |. Let us now rewrite B{t) = + B{t). From flT9|) and l\13b\} we 
expect that A — )■ as r — t- oo. We also have -B — )■ as r — t- oo. We shall now solve the 
system (I13atil36p through successive approximations. From the above assumptions we 
get for ffial : 

, 17 , 24 1 ~ (25)2 

A' + -A = Ye-^-(l -2B + ^ -...)- (23) 

-^e-''^l - 12B + ^^^^ -...)- - AB'\ 

We expect that A as well as B will be expressed as a series of e~'^"'^-like expressions for 
subsequent /x^'s so that they eventually all cancel each other. The first exponent to be 
cancelled is e~2'^ - therefore, neglecting higher order terms, we get 

A' + —A = — (24) 

4 5 ^ ^ 

with a solution 

A{r) = ^e-^^^ + ae-'^\ (25) 

where a comes as an integration constant for homogeneous equation. We can now 
rewrite fl23|) using A{t) = ||e~^'^ + ae~^'^ + A{t) to get: 

A- + EA = ?le-H-2B + <?4£ - . . .) - - 12B + <H|)! _ . . .) (26) 

4 5^2! ^ 5 ^ 2! ^ ^ ^ 

17 ~ 64 1 ~ 

-2ae--^B' e'^^B' - 2AB' - 45'1 

25 

Note, that B can not have a e~5'^ term because there would be no possibility to cancel 
^e~2^(— 25). Instead, this factor must combine to e"^"^ in order to cancel — |e~^'^. We 
can write 

e-H-fs - - = (27) 
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and solve it 



4 



Unfortunately, this is a wrong solution. One can see this by substituting it into fll3al) 
and trying to find next order of expansion - unbalanced exponents appear. This leads 
to a conclusion that A must also contribute to e~^'^ equation. It can be consistently 
done only by including e~^'^ term in A as well as e~t^ term in B: 



15, 



B{t) = \r + Ce-i' + 76"^^ + 5(r), 

A{t) = ^e-^ + + ae-x- + i(r), 

where C and D are yet unknown coefficients, B has been redefined and A has been 
defined for the next step of the method. After substitution into (11361) we get an equation 
for terms proportional to e~^'^ 

(25 D + 64 C + 16) = 0. (29) 

Equation f llSap provides us with yet another equation for e~^'^ term 

1 - 4^ = 0. (30) 

We therefore get C = D = — 1|. Everything is consistent and we can proceed to 
higher orders by guessing what terms should we include to cancel constant factors in 
equations for further exponents. It is interesting that we never arrive at any constraints 
for a or 7. This constants explicitly appear in higher coefficients of the expansion, but 
remain completely undetermined. 

Our heuristic reasoning is justified by theorem from Appendix A which ensures 
that the above formal expansion defines well behaved three parameter family of local 
solutions near infinity. The theorem requires that we interpret the integration constant 
E from ( 12T|) as a scaling constant and that we parameterize solutions with {0,7, E}. 
Below we present an expansion computed to one more order, translated into the original 
variable r = e"^: 

^, , 32 1 32 1 a 12 1 

A(r) ~ ^ -\ yf TT- + ..., (31) 

^ ^ 25 25 25 i?^ ' ^ ' 



2 



, , 1 , 11 7 3 1 

B{r) ~ - Ini? + -— + ^ + + . . . , 

^ ^ 4 4:Ri 8R^ 

where R{r) = Sr. 

We should remark here that the TN soliton exhibits exactly the same expansion 
with 7 = —j\/2, a = ||\/2 and S = ^^=. This feature enables us to define a mass of a 
WBH in terms of asymptotic integrals. 



5.3. Mass of the warped black hole 

The definition of a mass for the WBH is not straightforward because commonly known 
mass definitions (e.g. Komar mass [lOj) rely on asymptotic fiatness of a given solution. 
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Yet, our spacetime is not asymptotically fiat. It is locally asymptotically flat and we 
have to resort to other methods for mass computation. Abbott and Deser [11] have 
proposed a general definition of conserved charges. They replace asymptotically fiat 
background with a spacetime of an arbitrary asymptotics which is supposed to be the 
vacuum of the system. Next, they use its isometrics to generate conserved quantities. 
Below we present an application of this method to the case in hand. 

We choose the TN soliton as a background solution and decompose the metric to 
factor-out the background: 

Qap = QaP + hafi (32) 

One can define 

= i [g'^^H^^ + f^H"^ - g'^m^^ - g^^H"^) . (34) 

All indices are raised and lowered using the background metric, cja/B- The TN soliton is 
static, therefore it possesses a Killing vector field ^°'da = The charge associated 

with this isometry can be interpreted as a mass and equals to 

^ / [foV^ii"""' - K"-"VjU] dSi (35) 



where dSi is the spatial surface area element at large radius, the index t denotes the time 
coordinate index. Greek indices run through all spacetime directions while Latin ones 
run through space directions. The integral should be performed at spatial infinity. That 
explains why this method is sometimes referred to as method of asymptotic integrals. 
We can, therefore, use the expansion (13T]) for computation of the above integral to obtain 

4 ^ 5 16 

5 4 5 

This result is valid only if the scaling parameter E of the background TN soliton is the 
same as for WBH. Otherwise the integral is divergent. Note that S does not appear in 
the above formula. 

It is interesting that the WBH solution is parameterized by different number of 
parameters near th than near infinity. More precisely, near rn we have two parameters 
{f^yi^n} while near infinity - three {E,a,7}. Therefore there must be a constraint 
involving S, a and 7 which reduces the freedom at infinity. We determined this 
constraint numerically as follows. We chose th, shot the solution to infinity and fine 
tuned /3 to get S = 1 from the fit. Then we were able to fit a and 7. The results are 
presented in figure [Has a curve parameterized by th- Combining this knowledge with 
the equation (!36l) we can compute the mass of the WBH basing solely on th- 

One might wonder what will happen if we perform shooting from infinity towards 
horizon and take a and 7 lying outside the curve from Figure [H Using numerical 
simulations we looked at this case and observed that A > and becomes singular near 
r = 0. This is not a surprise, since the reasoning presented in the proof of existence 
of WBH (e.g. A can not cross 1 from below) is only valid when shooting forward, not 
backward. 
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Figure 1. The relation between asymptotic parameters a and 7 for WBH is plotted 
as a parametric curve, horizon radius rn being the parameter. It can be seen that for 
r/f — > the curve converges to {7, a\ of the soliton, 7 = — 1\/2 and a = f|\/2- 



5.4- Comparison to the 4+1 dimensional case 

In the 4+1 dimensional model described in pQ the black hole solutions are explicitly 
known. It is instructive to look at the parameter dependence between near horizon and 
far from horizon power series expansions and compare it to 6+1 dimensional case where 
only numerical relations are known. Let us begin by recalling the lower dimensional 
ansatz: 

ds^ = -Ae-^'dt" + A-^dr^ + -^r\^^{de^ + sin^ ed<P^) (37) 

+ ^r^e-^^{dij + cos ed(Pf. 

The A and B functions have an analogous interpretation as in 6+1 dimensional model. 
From the above ansatz we get field equations. Below we present their static version 

^(^5'yA + -r5'(4e-2^-e-8^) --(e-2^-e-«^) =0, (38a) 
3 3 

vA = -2A{l + r^B'^) + \ (4e-2^ - e-«^) , (386) 

5' = -2rB'\ (38c) 

where primes denote differentiation with respect to r. Next, we present KK solution 

5o = ^ln(l + £), ^0 = ^^^, e^'^ = A, (39) 

where 

r = 2m'/'p'/\l + ^)i/6 (40) 
m 

and m > is a free parameter which sets the scale (corresponding to / of ([1])). We also 
recall the black hole solution [T2l: 
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where 

A = p + 3M - VM2 + 2P^ S = p + M + VM2 + 2P^ (42) 

r = 42/3pl/3^1/6^1/2_ 

Here < P < 2Af . To see that this is indeed a black hole we can find local asymptotics 
similar to ffT2l) 



2(e-2/^ 



r 



4e-2.-e-«/ U-^)' 



A(r) ~ - (4e-2/5 - e-«^) ( — - l\ 
3 V^/^ / 



This time, though, we have explicit formulas connecting /3 and th with P and M 
^ 1 M + + 2P2 

^=3^^ 2 P ' ^ 

= 4^p2(M + VM2 + 2P2)(3M - VM2 + 2P2). 

Next, careful asymptotic analysis applied to fl38ati386|l gives us the following near- 
infinity expansions: 

N 16 1 a /8 16 2\ 1 

i^W~ilnP+^-(|a7 + 7^ + 0^..., 

where a, 7 are free parameters and R = 'Er represents scaling freedom. In contrast with 
the 6+1 dimensional case, free parameters appear just after the leading order term. Note 
also, that the next order depends explicitly on these parameters. For KK monopole we 
get a = — ^5 7 = i' ^ — 2m' black hole we have 



8 M + VM2 + 2P2 _ M - VM2 + 2P2 1 

Therefore P is simply the scale factor. It is evident that two parameters {M, P} are 
translatable to {(3, r^} on one hand and to {a, 7, E} on the other hand. 

To compare these results with the 6+1 dimensional case we set P = ^ and get 
S = 1. One can now see that th depends solely on M > | and that for M = | we have 
th = - black hole degenerates to the KK monopole. (3 is also fixed by the choice of 
M so we can parametrize the black hole up to scaling by rj^ alone. 

Next, elimination of M from (H6|l gives us a surprisingly simple relation between a 
and 7 

7 = -^-- (47) 

9 a 

It can be easily seen from Figure [1] that in the 6+1 dimensional case the corresponding 
relation is very different. 

In the end we present mass of the 4+1 dimensional warped black hole computed 
using the same method as for 6+1 dimensional case: 

M4+i = -^(8 + 327 + 9a). (48) 
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After substituting (H^ . we get 
6. Numerical results 

In order to investigate nonlinear perturbations we have resorted to numerical 
simulations. Using second order finite difference code we have solved the equations 
f l4all4dj) for several families of regular initial data representing perturbations of Taub- 
NUT soliton ([1]). The overall picture does not depend on the specific choice of a family. 
The results shown below were generated for initial data of the form (using momentum 
variable P = e^A-'^B): 

B{0,r)=Bo{p{r)),PiO,r)=p(^^^ e^. (50) 

The scaling parameter of the initial soliton has been fixed to = |. Parameters a i tq 
are fixed for given examples, whereas the amplitude p varies. 



6.1. Small perturbations 

It has been observed that for small perturbations, by which we mean small value of 
the control parameter p, the perturbation is being radiated to infinity and the solution 
settles down to the Taub-NUT soliton with the same parameter / as the perturbed one. 



One might wonder why / is not altered by the evolution. If we recall that / is 
connected with the scaling constant by relation T, = we see that changes in / would 
change scaling. The mass formula fl36|) becomes then divergent - it would require an 
infinite energy to change scaling, whereas our perturbation is of finite energy. 

An interesting aspect of the relaxation process is connected with the absence of 
term in expansion (13T|) of A. If one examines (l¥aj) closely it can be seen that the 
compact-support perturbation imposed on the momentum variable instantly introduces 
r""^ term in the expansion of A for large r. If the perturbed soliton is really stable, this 
term should eventually vanish. More precisely, from ( l3T|) we know that the function 

d{t,r) = r^{A{t,r) - Aoir)) (51) 

represents the coefficient in front of term. Figure H] presents the behaviour of this 
function. 



6.2. Late time tail 

In order to better understand the relaxation process we transformed the linear equation 
([9]) into single mode wave equation 

- -r^u{x) + V{x)u{x) = X\{x) (52) 
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Figure 2. Asymptotic stability of the Taub-NUT soliton in 6+1 spacetime dimensions. 
For initial data (|50|) with a small amplitude {p = 0.02, — ?>, a = 1) we plot a series 
of snapshots of the function B(t,r) (where t is central proper time). The dashed line 
shows the unperturbed TN soliton. During the evolution the excess energy of the 
perturbation is clearly seen to be radiated away to infinity and the solution returns to 
equilibrium. 



using transformation 



V^p" (p + 4)= 



The potential has the following asymptotic behaviour 

V{x) = -, + ^-^ + -,+0{^). (55) 

Oil Ob OC tC 

According to work of Ching et al. [13] and Bizoh et al. [H] the solution should exhibit 
asymptotic power-law tail 

(56) 

We are aware of the possible existence of nonlinear tail, however our numerical 
simulations are in good agreement with ( 156|) . Furthermore, the amplitude of the tail 
scales linearly with initial data, which suggests that even if nonlinear tail is present it 
decays faster than the linear one. 

Our simulations strongly suggest that the TN soliton is nonlinearly stable against 
small perturbations. However, just like its lower dimensional cousin, it is being 
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5 10 15 20 25 30 35 40 

t 

Figure 3. The convergence to the TN sohton. For the same initial data as in Figure 
[5] we plot ln|P(ro,t)| as a function of time for tq — 2.0. Asymptotic power-law tail 
can be seen. 




5 10 15 20 25 



Figure 4. Relaxation of the term in asymptotic expansion. For the same initial 
data as in Figure [5] we plot series of snapshots for d{t, r) as defined in ([?T|) (t is central 
proper time). The dashed line indicates zero. During the evolution the coefficient 
goes to zero behind a wave travelling towards infinity. 
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destabilized by sufficiently strong perturbations. Warped black holes are important 
in this process. 

6.3. Large perturbations 

We have found that for sufficiently strong perturbations the inner part of the soliton 
collapses and horizon is being formed for r = r^^ > 0. The external part of the solution 
settles down to a WBH. The scaling parameter / is preserved. 




Figure 5. Instability of the Taub-NUT soliton for large perturbations. For initial 
data ([50)1 with large amplitude [p — 0.1, rp = 3, cr = 1) we plot a series of snapshots 
of the function A(t,r) (where t is central proper time). During the evolution A(t,r) 
drops to zero near r — rn — 2.174 which signals the formation of a horizon there. 
Outside the horizon the solutions settles down to the warped black hole (dashed line) 
with M = 0.612. 



7. Conclusions 

W have shown that the TN soliton in 6+1 dimensions is classically stable just like the 
KK monopole. Therefore the stability is neither dimension specific nor supersymmetry 
dependent feature and its origin should be sought elsewhere. We have also shown that 
for sufficiently strong perturbations a collapse to the warped black hole occurs. The 
black hole properties have been studied in detail and the proof of global existence has 
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been presented. Still, there are open questions concerning the model. The first one 
is connected with the critical behaviour on the threshold of black hole formation. In 
4+1 case the discrete type II self similarity was observed and it might be interesting 
to compare those results with the observations concerning the case in hand. The other 
interesting feature that can be examined is the appearance of quasinormal modes in 
asymptotic pointwise convergence presented in the Figure [31 
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Appendix A. 

This appendix is devoted to the formal proof of the existence theorems for static 
equations. Let us begin with the theorem which has been proven in [9] (Proposition 
12): 

Theorem 1. Suppose a system of differential equations for m+n functions u = 
(mi, Um) and v = (f i, f„), 

t^ = t'''fi{t,u,v),t^ = -\iVi + t'''giit,u,v), (A.l) 
dt dt 

with constants Aj with Re{Xi) > and integers fii, Vi > 1 and let C be an open subset 

of such that the functions f and g are analytic in a neighbourhood of t = 0, u = c, 

V = for all c G CC. Then there exists a unique m-parameter family of solutions of the 

system liA.l]) with boundary conditions {u,v) = (c, 0) at t = 0, defined for c E €. and 

\t\ < to{c) with some to(c) > 0. These solutions satisfy 

Ui{t) = Ci + 0{t^^),Vi{t) = 0{t''^), 

and are analytic in t and c. 

Now, we cast static equations into the form: 

{Ar^'B')' = -Ar^AB'^ + -r^ (e'^^ - e'^^^) , (A.2) 

A = ^{1 + r^B") + - (24e-2S - Ae-'^"") , (A.3) 
S' = -ArB'^. (A.4) 
The first two equations are independent of the third. This allows us to consider only 
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Appendix A.l. Existence near r = 

We will try to prove the relations ([7]). In order to do this, we perform a transformation: 
1 - AB' I- A 

According to our supposition on local behaviour we should get: 

uJi{r) ~ h,uj2{r) ~ 2h,ujz{r) ~ 86^. 
Now, equations (]A.2|A.3P take the form: 

ruj[ = -2uJi +UJ2 + r^/i, 

= I2uji - 6u2 + r^f2, (A.5) 

rco'^ = 4 {ul + 12^1^ _ 2^3) + r'fs, 
where fi are polynomials in w, ^,r^. We need to perform another transformation: 

Ul=Ui- Vi, 

U2 = 2ui + 6vi, (A. 6) 



1^3 = ^2 + 8{3vl + 
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to finally get 



rv[ = —8vi + r^/ii, (A. 7) 

rv2 = -8v2 + r^h2, 

where g, h are polynomials in u,v,^,r^. This system has the form ( lA.ip . thus from 
Theorem [1] we acquire 

Theorem 2. There exists an unique, one-parameter family of regular solutions for 
\A.!^) near r = analytic in b and r such that 

B{r) = br^ + O(r^), A{r) = 1 - 86V^ + O(r^), 

(A.8) 

Appendix A. 2. Existence near th '■ A{rH) = 



We need a transformation which casts ( 1A.2tlA.3l) to the form ( lA.ll) . around th > ^ '■ 
A{rH) = 0: 

1 - 

p = r -rH,Ui{p) = r,U2{p) = , (A.9) 

A AB' 

OJi = — , UJ2 



We get 



P P 



pu[ = p, 

p4 = P— , (A.IO) 



pu[ = -ui + Fi{p) + pfi, 

pu'2 = -UJ2 + F2{p) +P/2, 



where Qi are polynomials in ui,U2,-^,^- Those polynomials are regular around 
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where /i and /2 are polynomials in wi, a;2, and 

i^i(p) = ^ (24(1 - 2u2) - 4(1 - 2u2f) , (A.ll) 

^2(p) = ^(l-2M2-(l-2n2)'). 

Next, transformations 

uji = Vi + Fi (A.12) 

give : 

pVi = -Vi + pgi{u,v), (A. 13) 

j_ j_ 

Ui = rH,U2 = f3,Vi = 0. According to Theorem [1] we get: 

Theorem 3. There exists an unique two-parameter family of regular solutions of the 
system /lA.^A.!^) around th '■ A^rn) = 0, rn > 0, analytic in (3, th and r, so that: 

B{r — th) = P + 0{r — rn), A{r — th) = 0{r — th)- 
Appendix A. 3. Existence near infinity 

Let us consider a coordinate change u = We then take p = T,u, S > 0: 

^(P) = |(P^-P^^) + 21(P)P^ 
5(p) = _lnp+lpio + ®(p)pi^ 

S{p) = ^. 
u 

Some more transformations are required: 

g-ipiO(i+4p^'8(p)) ^ 1 _ 1^10^1 ^ 4p5a;(p)), 

-(^)-|- 

Those transformations leave us with the set of equations: 

p— = pf{Ql,u,S, Z,p), 
dp 

dS 

p-f- = P 0, 
dp 

P^ = P^(P), 
dZ 

p— = -16Z(p) + pg{Ql, u, S, Z, p), 
dp 

where /, g are regular for a; — )■ 0, 21 — )■ a, w — )■ 7, — )■ S > 0, p — )■ 0. Those equations 
are of the form (lA.ip . We now get from Theorem [TJ 
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Theorem 4. There exists an unique three-parameter family of regular solutions for the 
system around infinity with asympotics B{r) ~ Inr, analytical m S > 0, a 

and 7, so that: 

4/ N 32 1 32 1 a ^, 1 , 



1 1 1 7 1 

B(r) = -lnR+ + + 0(—r' 



where R = T,r. 



Appendix B. 

We present an explicit form of V{p) introduced in the proof of linear stability : 

^(P) = 7^7?7 nnA^ 9 \ r .,9X9/ Jn >< fs^V + 60/^" + 502/V° + 

+ 2250/V^ + 5665/V^ + 7500/ V + SGp^P^Q^ - Sp^PQ^ + AOOOpH^ - ^Ip'PQH + 
+ 504p^P^Q^/ - mp^PQH'^ + 2844p^p6Q^/2 _ ggOp^PQ^p ^ %2mp^P'^QH'^ + 

+ 13140p^P^g^/2 - 1095p2pQ^/^ - 900pPg^/^ + lOSOOpP^Q^/^ - 'i^^PQH^ + 3600P^g^/^) 

where 

/ /^(p + 4o y 
V(p + 20V ' 

Q=(/V(p + 2/)2(p + 4/))^. 
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